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ABSTRACT 


A simplified form of the three-degrees-of-freedom nonlinear 
equations of the 3000 ton Surface Effect Ship (SES) was lin- 
earized around an equilibrium point when a thruster failure 
or thruster reversal occured. Optimal control feedback gains 
were calculated with the linearized system and applied to the 
nonlinear system. Evaluation of this design was based upon 
both comparison of the real-time human operator response and 
additionally a state feedback design based upon a trial and 


error procedure. 
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ieee PRODUCTION 


Over the past years a lot of attention has been given to 
the Surface Effect Ship (SES) due to its advantages over con- 
ventional ship designs [Refs. 1, 2, and 3] 

Recent investigations regarding water jet propulsion SES 
have been conducted on the computer assisted mode of operation 
either in a normal mode of operation or during a variety of 
failures [Ref. 3]. Failures considered include, among others, 
complete loss of thrust and reversal of thrust of a nozzle. 

The purpose of this thesis was to investigate the possi- 
pelaty of optimal maneuvering of a 3000 ton Surface Effect 
Ship during the above mentioned emergencies. 

Ms discussed in [LRef. 3], turn maneuvers for the 3KSES 
are the result of turning moments generated by modulating the 
waterjet thrust vectors. This modulation is achieved by a 
combination of thrust amplitude modulation (controlled by 
varying the flow rates of the four waterjets to create a 
net differential thrust, i.e., by varying the power turbine 
speeds of the individual propulsion engines), and thrust 
vectoring (controlled by varying the directions of the four 
waterjet deflection sleeves). 

In the computer assist mode, the relative position of the 
helm determines the angular rate at which the turn is nego- 


tiated. Full helm position corresponds to the maximum safe 
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turning rate for a given ship speed. Calibration of the helm 
input is done within the control computer, which then gen- 
erates commands to the individual control effectors to produce 
the commanded turn rate. 

In the manual mode, the relative helm position generates 
commands for the direction nozzles directly. A common position 
command drives all four deflection sleeves. Because of that, 
the rate of turn which corresponds to a particular helm posi- 
tion is dependent on the hydrodynamic yawing moment which is 
a function of ship speed. The operator can also generate a 
thrust differential turning moment in the manual mode by 
proper selection of the propulsion engine control levers. 

To relieve the operator of the work load related to moni- 
toring and maintaining ship heading over extended periods of 
time, a heading control autopilot is included in the design 
of the heading control system. From an operational point of 
view, the heading autopilot is limited to a "heading hold" 
function, and is used only in conjunction with the computer 
assist mode of operation. The net result of this operation 
is that the heading autopilot replaces the operator's hand on 
the helm, while the helm is still the primary control input 
deyice for the control of ship heading. 

This is accomplished by direct feedback of the yaw rate 
to a processor which also gets an error channel from the auto- 
mieeet through a control filter using proportional plus inte- 


gral control of the heading angle. 
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The fundamental optimality requirement used in [Ref. 2] 
is concerned with minimizing fuel and effector duty cycle 
while executing prescribed maneuvers. This was done by 
Minimization of a function of the state variables. 

The approach used in this work was to achieve results by 
Minimizing a functional of the form ptt 


ceo 


to compute the gains for the state feedback control (see 


(x Qx + wR) dt 


Appendix A). 

The resulting linear regulator design was tested on the 
nonlinear system to compare the performance of the craft 
during a failure mode with the performance of a human operator 
using a real time simulation [Ref. 4]. The performance of the 
design during a failure was also compared with a corrective 
maesorithm for control found by simple trial and error pro- 


cedures. 
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fener or MOlLON OF THE THREE-DEGREE- 
OEP REeDOM CSbOr) MODEL 


The three degrees of freedom to be used in this analysis 
are: surge, sway and yaw [Ref. 5]. Therefore, the following 
assumptions and approximations were made: 

i Roll is zero, 

2. Pitch and draft remain constant. 

3. The origin of the ship's coordinate system is located 
at the center of gravity of the ship. 

4. At time = 0, the ship's coordinate system and the 
fixed (geographical) coordinate system have the same origin 
and the same coordinate axes, i.e., 


eee og Eo Yay 
fe 0 (see Fig. 1). 


O 
5. At time = 0 the ship is heading in the positive Xray 
Saeeection. 
6. Drift angle subtraction effect on thrust angle is zero. 
Then, aS a good approximation, the equations of motion 
are: 
Surge: m(u-vr) = X Gi) 
Sway: m(vtur) = Y (2) 
Yaw: T_D=N C3) 
Navigation: 
XS = ucost-vsinw (4) 
Yo = uSinwt+vcosy (5) 
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Fig. 1 - Definition of Coordinate System 


NOTE: Rudder condition for right turn is shown. All runs 


start at geographic origin. 
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oD 


2 -l,-vV 
Cea an c? Ch) 


mathe above equations: 


m= mass of the rigid ship 


I. = moment of inertia about the z axis 
u = V .cosB = velocity in the x direction (Surge) 
Vv = V sing = velocity in the y direction (Sway) 


r= w = angular velocity about the z axis 
vw = heading (yaw) angle, Fig. 1 
B 


=Gritt angle, Fig. 1 


6 = rudder angle, +é produces a left turn 
Also: 

a fx dt Ee fy dt 

us = fudt Vie =a dt 

a = frdt Y= ava te 


The right-hand side of the equations of motion have been 
represented by X, Y, N. These symbols represent: summation 
Merorces (X, Y) and moments (N). Each summation consists of 
many terms. For the initial simplified model, only the few 
terms that are major contributions have been retained. Those 
terms are considered sufficient to produce a simulation which 
would be in close agreement with obtained data. At a later 
time, the model may be improved by adding more terms and ad- 


justing coefficients. 
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ee eee dem omeotethe rightrhand=sidées, these force and 


moment summations become: 
2 


X = GTcosd-Ch yu (8) 
Y = GTsiné-Chv|v| ce) 
he = 


-(G6Tsiné)1 +C,,, 


viv|4 (10) 
W 


where (see Fig. 2): 
eue= gross thrust in lbs. 
5 9) ees = angle of all thrust vectors (starboard and port) 
in radians 


= coefficient of drag in the direction of surge 


DX 
Chy = coefficient of drag in the direction of sway 
Lo = thrust level arm in feet 


1. = sway-drag level arm in feet 


Also, note that the drag forces have been chosen as lumped 
parameters acting at the two points P & P', i.e. P & P' are 
equivalent force centroids. 

The length Lo is the distance from the center of gravity to 
the stern, and 1. is the distance (in the x-direction) from the 
center of gravity to the points P & P!, 1) is taken as a fixed 
Parameter (it could be charged by various conditions of ship 


Weeding) but 1. is adjustable in order to be used as a para- 


meter to make the model fit test data. 

The coefficient Chy in eq. 8 is a lumped parameter equiva- 
tent for the hull friction effect (slender body theory). The 
term ate models the hull drag, aerodynamic and hydrodynamic 
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effects. The bubble drag is not included in the initial 
model. 


In like manner, the coefficient C in eqs. (9) and (10) 


IBS¢ 
models the cross flow drag effects. Note also that the added 
mass term has been temporarily omitted and will be added 
later, especially in the 3KSES study since its contribution to 
the XR-3 model is almost nonexistent [Ref. 6]. 

SMstitution of equations (8), (93) and (10) into the 


equations (1), (2) and (3) provides the basic equations of 


motion for the three-degrees-of-freedom model: 


i = + (GTcosé - C ae) + vr (11) 
Dx 

Vv = = (GTsiné - Coyv ivy? - ur Gi?) 

mee 1 : 

r = a (-(GTsing)1, + Cpyvi|v] 1.) (13) 


Bawations (11), (12) and (13) are applicable to any SES 
craft. Only the particular parameter values to be used in 
the equations would depend on the particular craft (XR-3 or 
Sob.) . 

In the case of the 3KSES, the equations take the form 


Maee Fig. 3): 


=) 


a2 (-C, u°+(T74+T8+T9+T10)coss ) + vr 


= 1 , 
a ((T7+T8+T9+T10)siné -C,v|v[)-ur 
(Cry (ww) v |v | + [T7coséS,+T8cosés, 


- T3coss 8, - T10coséS, 


Peles ing ) COO) ] + Ay suv (ww) ) 


="(T7sindtresind + T9sind 


lke 





7 
cry 


era & 


DY 


viv 


dS, 








(s~ 





per 


ucosw-vsinw 


Xe 
It 


usin) + vcoswt 


Ke 
O 
tt 


For the purpose of this study the thruster's angle 6 was 
the same for all thrusters except the starboard one, which 
in the case of thrust failure was considered 0° (thruster 
oriented straight ahead) and in the case of thrust reversal 
it was considered as being 180° (thruster oriented straight 
astern). 

In the first case, the thrust of the failed thruster was 
Set to zero pounds and in the second case the thrust was set 
to full value (98000 pounds for 56 Kts operation speed). That 
means that either in the case of thrust failure or thrust 
reversal of the starboard thruster, the rest of the thrusters 
are ganged together and act in the same direction. 

Therefore, the angle 6 can be redefined as 27, 28, Z93 and 
Z10 as shown in Fig. 3 to correspond to the appropriate thrus- 
ters and in both cases Z7=Z82Z9 with 210 being either 0 or T 
radians. 

In order to linearize the equations, the following assump- 
tions were made: 

1. The maximum rudder angle is less than or equal to 30 
degrees and for most high speed manueuvers the rudder angle 
Stays below 15 degrees. Thus, by letting cosé*1 and sind*6 
(where 6 is in radians), the difference from the actual values 


memo = .2618 rads = 15° is 
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a. cosl5° = .966, difference = aaa x 100%=3,5% 


~2618-,2588 
2.010 3 


De sans? .2588, difference = 


roo = 1.154 

Differences with these magnitudes should not introduce 
large deviations from the linear design when applied to the 
nonlinear equations. Thus in the analysis eee was restricted 
to 15 degrees or if a exceeds 15 degrees it probably will 
be of short duration and will then not introduce large errors 
in the design values for control. 

2. The yaw angle » was considered to be less than 10 
degrees which agrees very well with the results obtained from 
Seheesimulation of the craft. 

3. Even though the sign of the sway velocity changes de- 
femaenes On the particular situation occuring, the term viv! 
was considered as ve 


With the above assumptions, the equations for the 3KSES 


become: 
f, = = = (-Cpu +T7+T8+T9+T10) +tvr 
el: 2 

fe = ((T7+T84+T94+T10)6—-C...v) - ur 

2 m DY 
io-  ~ = (ec vy (ww) 4+T7S #+T8S.-T9S 

3 1, iba 3 IL 

~ TLOS,,-(T7+T8+T9+T10)5 (OO) +A, ,uv (ww) ) 

me Y= 


Za: 





; = TO ae ie 
a y 


Defining the state variables 


xX, Fu-u 
PTY Ve 
A 
X_=P-V 
A 
ce! Yo, 
A 
Xp =X a 
- 
Bo) -0o "coo 
: 7 
Wj et es 
U5227-27 | 
Then, xy =u 
Se 
X 4 = 7 
Xi, = ¥ 
Xe = a 
ists : Yo 


The linear state equation coefficients are: 
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r = 7 
men 2 
O 
V = V 
V=V : 
O 
Bay : of, Po 
ox | oy 
L 
m 
0 
36 es Ese. 
r=r ‘ 
O 
-2Cpyv : ~2Cp vv 
m. - 7 m 
V=V 
O 
—4)5) = TUR 
u=U 
O 
ae = aa = (0) 
ox dy 
27, nih 6 
e 27227 = 
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In standard form the system becomes: 


~“ = Ax + Bek where 


Zo 
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by linearizing around the equilibrium point Cua Vo: 


ee yao the equations become: 
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The above equations were used to determine by computer 
Simulation the required gains for the state feedbacks using 
the linear regulator solution [Refs. 7, 8 and 9]. 

In order to obtain the equilibrium point, the nonlinear 
equations were used by setting u=0, v=0, r=0 and determining, 
after a failure had occured and the system being corrected as 
far as was possible by the trial and error algorithm (see 
Chapter III), the "steady state" values for the above men- 
tioned variables. 

Two "steady state" runs were made, one with the starboard 
thruster failed and the other with the starboard thruster 
reversed. 

The "steady state" values obtained are listed in Tables 


imeand Li. 


TABLE 1 


Steady State Values (Thruster Failure) 


Surge Velocity Sway Velocity Angular Velocity Yaw Angle Track 


Deveson 
feet /s) CE) sD (Caady s)) (rad) Cp) 
7 ay, 1.86 = pinata =a 54 ile Ome ss: 


Bert Lhrust 
pT) 


87600 
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TABEE IL 
Steady State Values (Thruster Reversal) 


Surge Velocity Sway Velocity Angular Velocity Yaw Angle Track 


Deviation 
jm c/s) Cie/ Ss) (rad/s) (rad) Chie) 
55.66 Pee? aena7sa07 > = 07 09 Par 


For this study, the steady state value achieved by the 
Simulation with the thrust failure were used for the develop- 
ment of the linear system and the subsequent use of the opti- 
mal control scheme. The choice was quite arbitrary but, 
nevertheless, the results show that the system does behave 
well, no matter which one of the two emergencies (thrust 


Mmeiir~e or thrust reversal) occurs. 
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eee een ei en ROR DESIGN 


A. INTRODUCTION 

Taking into account the impressive speed of the surface 
Seeecrt Ship, it is evident that there is a need for controlling 
it either during straight runs or during an emergency which 
would be manifested by the loss of thrust of one thruster or 
even by the thrust reversal. 

As was seen from the results of computer runs with the 
model, the behavior of the ship during a failure was pre- 
dictable but nevertheless disturbing because the ship veered 
quite rapidly from the track that it was following before 
the failure. 

Therefore, the goal of the control procedure was twofold. 
First, it was required to maintain control of the ship and 


second to maintain that control in some acceptable manner. 


B. ANALYSIS 

As was evidenced (see Fig. 4), when the outer starboard 
thruster experienced a complete thruster failure, the ship, 
without any corrective action with respect to nozzle movement 
or thrust variation of the remaining jet nozzles, veered to 
the right and actually entered a turn. 

By analyzing this behavior, it was seen that the para- 
meters changing the most were the yaw angle w and the deflec- 


tion from the original track which was assumed straight ahead. 


od 





Intuitively then, it was evident that from the available 
Variables, w and Ve would be most eligible to control in order 
to bring the ship back or close to the original track. Be- 
cause it was quite clear that, if the heading (yaw angle) was 
Minimized as well as the deviation from the track, the ship 
would follow again the original track or a track that was 
Semaiilel to that one. 

Therefore, the corrective action chosen was nozzle de- 
flection as a function of the yaw angle and the deviation 
from the original track. 

Furthermore, since positive nozzle angle drives the ship 
to the left, by selecting deviation to the right as being 
positive and yaw angle to the right being positive to the 
right, it was quite clear that the nozzle action would only 
involve yw and vie without any sign change. 

Thus, the problem reduced to the point of determining the 
appropriate gains for w and Y, and the rudder action took the 


= Le 
ZOrm Z KY ROY: 


feet Al AND ERROR PROCEDURE 


As a beginning, K, was arbitrarily selected to be 1000 


I 


and K being 5, It was seen that the action of the rudder 


was too heavy and the movement of the ship was quite erratic 
and at the end became unbounded. 
Therefore, both gains were decreased to K,=8 and K,=.035, 


which gave a more acceptable performance (Fig. 14). 
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The thus established rudder algorithm proved to be adequate 
even when a reversal failure of the starboard thruster was 
mmeiacted (Fig. 21). 

As a further improvement, the idea of modulating the 
thrust of the left thruster during a failure was brought into 
effect. That is, the thrust was varied to smaller values by 
the amount of the deflection of the ship off the original 
mpack (Fig. 20). 

The above scheme was used for all the runs presented in 
Mees. 14 through 27. 

In order to evaluate all the above corrective actions, 
the performance indices described in Chapter IV were used. 

As can be seen, the ship actually moves down the track 
with a large drift angle. That is, although a corrective 
rudder action is maintained, the ship does not turn into a 
circle and also it does not veer to the left or to the right 


although there exists a sway velocity and asymmetrical thrust. 
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IV. LINEAR REGULATOR DESIGN 


fe LNTRODUCTION 

The theory behind the linear regulator is well established 
and need not be repeated here. For a brief discussion, see 
Appendix A, 

After the linearization of the system's equations around 
an equilibrium point, the IBM/360 was used to obtain the 
solution of the matrix Riccati equation, using the method 


presented in [Refs. 8 and 9]. 


B. ANALYSIS 

It was evident that in order to achieve a parallel track 
for the craft after a failure had occured, the deflection 
from the original track should come into play. Also, in 
order for the craft to maintain this parallel track, the yaw 
angle should be minimized as well as the angular velocity 
m@pout the Z axis. 

As with the trial and error procedure the corrective 
action chosen was ganged deflection of the operating nozzles 
Smemthnrust control of the port thruster to compensate for the 


fmabbed Starboard thruster. 


Seer oREORMANCE INDICES 
One of the basic problems of modern control theory is 


the translation of system specifications, often in such 
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Subjective terms as "good rise time with reasonable overshoot", 
into the form of a performance index. Furthermore, after the 
selection of a performance index, the designer has to make 
an initial guess for the weighting of the state variables and 
adjust them after comparing time histories of the controlled 
system. 

For the purpose of this study, the following performance 
index was selected: 


w= fee (x Qx al 


ene) aclets 
O ~o lu ON 


Concerning the weighting of the variables, considerable ex- 
perimentation was done. In the final form of the above 


performance index the weighting matrices Q and R were: 


0 60 0 0 0 0 
Diy ve 0 0 0 0 
0 oO 20000 oO 0 0 
Q = loo 0 Sane 0 
0 0 0 0 0 0 
0 0 0 a @7 {eG 
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The final weights used reflect the importance given to the 
yaw angle, the angular velocity about the z axis and the 
craft deviation from the original track. 


The optimal control inputs obtained are: 


0 0 SOCOM A ialeley, 0 0 


— 
TT 
ax 


-.0316 .684 sheen) Ee Oe 0 S66 


The response of the nonlinear system with the above 
feedback is shown in Figs. 28 through 41. 

in order to compare the performance of the craft under 
the trial and error algorithm, the optimal algorithm and the 
real time simulation, the following two additional perfor- 
mance indices were used: 


J, = sf y. dt and J. = stf (Z7) “dt. 
O O 
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V COME oon Or RESULTS 


A. THRUST FAILURE MODE (OPTIMAL ALGORITHM VS. TRIAL 
AND ERROR ALGORITHM) 


Observation of the sea track (x. VS. Se CF he oiadliiineona | 2.0 ) 
indicates that the trial and error algorithm is much less of 
the desired control algorithm when compared to the optimal, 
because it appears that the trial and error algorithm does 
not react fast enough to stabilize the vehicle whereas the 
optimal does so immediately. The trial and error algorithm 
does seek out the steady state value of the track deflection 
but the comparison of the two tracks is dramatic in showing 
how much better off is the optimal algorithm in controlling 
the craft during a thrust failure. It must be kept in 
mind in making this observation that the trial and error 
design used in this comparison was not developed to its full 
capabilities. As stated earlier, the trial and error design 
was intended only as a means of establishing the linear 
system design. 

From the time history shown in Figs. 14 and 28, the vehicle 
for the same failure scenario reaches a steady state point with 
the optimal control in 22 seconds whereas with the trial and 
error control the steady state point would be reached some- 
where around 200 seconds or even greater, due to the lightly 


damped control system. 


Sid) 





As can be seen from Figs. 20 and 35 the port thruster is 
much more active with the trial and error algorithm. The 
range of modulation covers about 14000 pounds whereas utiliza- 
tion of the optimal algorithm gives a total thrust modulation 
of 2 pounds. This does not reflect exactly what happens to 
the port thruster, because as stated earlier in Chapter II as 
soon as the failure occurs, the port thrust drops to 90% of 
its full value. 

It should also be noted that the trial and error solution 
seeks the steady state value of the port thruster obtained by 
the optimal solution. 

Pemparing Figs. 18 and 32 it can be seen in Fig. 32 that 
the forward and side thrust as well as the yaw moment reach 
a steady state within 10 seconds after the initiation of the 
thruster failure when the optimal solution is used. Instead, 
when the trial and error solution is used, the craft as shown 
in Fig. 18 seeks the same steady state but even at 100 sec 
after the initiation of the failure it has not stabilized. On 
the average, the values obtained from both algorithms are 
approximately the same. 

However, it is readily apparent that the trial and error 
control system has poles close to the jw-axis and as a result 
exhibits a response of a lightly damped system whereas the 
optimal control system responds with much faster time constants, 
Which would indicate that this system has poles at a greater 


distance from the jw-axis. 
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Pee loRUSLT REVERSAL MODE (OPTIMAL ALGORITHM VS, TRIAL 
AND ERROR ALGORITHM) 


All the previous comments that were made for the thrust 
failure mode also apply to the case of a thrust reversal for 
the starboard thruster. However, due to the increased 
severity of the forces acting during the thrust reversal and the 
greater magnitude of the yaw moment (see Figs. 32 and 39), 
the optimal control exhibits a slightly longer settling time 
Mates. 28 and 35). 

During the thrust reversal scenario, the absolute magni- 
tudes of the state variables is larger than those experienced 
during the thrust failure case. 

Comparison between the optimal solution time histories 
and trial and error solution time histories again confirms 
the fact that optimal state variable feedback provides a more 
highly damped system, with steady state values being reached 
within 20 seconds from the initiation of the failure (see 
fees. 36, 38, 39 and 40). The trial and error solution is 
shown in Figs. 22, 24, 25 and 26 which again exhibits a 
lightly damped time history which would not be acceptable 
either from a mechanical or ship's efficiency viewpoint. 

It is observed (Fig. 23) that the trial and error solution 
requires the operational thrusters to be extended to their 
mechanical limits of +.5 radians even at the time =100 seconds 
Moant, whereas the optimal solution in Fig. 37 has a peak 
value for the thrusters of approximately +.4 radians and 


never requires them to extend to their mechanical limit. 
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Fig. 4 - Sea Track, Thruster Failure, No Correction 
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Fig. 9 - Sea Track, Thruster Reversal, No Correction 
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Fig. 14 - Sea Track, Thruster Failure, Trial and Error Solution 
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Fig. 15 ~ Surge and Sway Velocities Vs. Time, Thruster 
Failure, Trial and Error Solution 
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Fig. 18 - Forward Force, - Foree and Turn Moment Vs. Time, 
Thruster Failure, Trial and Error Solution 
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Fig. 21 - Sea Track, Thruster Reversal, Trial and Error 
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Pig. 22 - Surge and Sway Velocities Vs. Time, Thruster 
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Fig. 24 - Yaw Rate, Heading and Slip Angle Vs. Time, Thruster 
Reversal, Trial and Error Solution 
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Fig. 26 - Surge, Sway and Yaw Accelerations Vs. Time, Thruster 
Reversal, Trial and Error Solution 
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Fig. 31 - Yaw Rate, Heading and Slip Angle Vs. Time, Thruster 
Failure, Optimal Solution 
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Fig. 35 - Sea Track, Thruster Reversal, Spend soluclLon 
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Fig. 38 - Yaw Rate, Heading and Slip Angle Vs. Time, Thruster 
Reversal, Optimal Solution 
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The sea track (see Fig. 35) of the optimal solution 
becomes stabilized at approximately 18 feet to the right of 
the original track (due north), with a slight tendency to 
Mere tO the right. This is not of significant consideration 
due to the relatively small magnitude of drift and the fact 
that during the conduct of an actual emergency situation 
proper corrective action would be also taken by the crew of 
the ship in the first 10 to 15 seconds of the emergency. 

It 1s important to keep in mind that the ship's response 
for both the trial and error and the optimal solutions was 
the result of totally automatic operation and any human 
operator inputs were not considered. 

The deviation to the right obtained by the use of the trial 
and error algorithm at 100 seconds as shown in Fig. 21 has an 
average value of about 18 feet and is fluctuating around this 
value by a magnitude of 6 feet which would be more than enough 
mo discrupt normal ship activity. 

In addition, the situation 1S more aggravated when surge, 
sway and yaw accelerations (see Figs. 26 and 33) are taken 
into account because even though they are small, they are 


Serie sienificant after 100 seconds. 
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MmemocolurARISON OF THE OPTIMAL SOLUTION AND THE 
feo 1iME HUMAN OPERATOR RESPONSE 


As stated in the introduction, a real time simulation 
of five-degree-of-freedom Surface Effect Ship was available 
[Ref. 4], utilizing a human operator to input only ganged 
thruster action, the thrust being constant. 

Comparison of the optimal solution and the above men- 
tioned simulation gave the following results during a thrust 
reversal scenario: 

1. Initial recognition of the emergency by the human 
operator did not occur until approximately three seconds into 
the emergency, aS .compared to the optimal solution with an 
instantaneous recognition of the failure. Consequently, 
as shown in Table IV, the absolute values of all the outputs 
reached by a human operator controlled SES are much greater 
than those observed in the optimal control, given the same 


Maitial conditions. 


siya ay ee 


Comparison of Optimal Solution and 
Human Operator Response 


HUMAN OPERATOR OPTIMAL SOLUTION 


Time for first overshoot(sec) 39 Le 
Deflection from original 100 Le 
frack (ft) (worst case) 

Time to settle back(sec) 78 Sis 
Binal deflection(ft) 0 We 
Turn rate(rad/sec) 007 BOG 


19 





Specific examples are as follows: 

1. Maximum track deviation obtained by the ship with the 
optimal control in effect was approximately 18 feet, whereas 
with the human operator maximum track deviation varied from 
mo teet to 100 feet, depending upon operator response and 
familiarization with the failure as well as the ability of the 
Operator to recognize the onset of an emergency situation 
through the visual cues provided by the real time simulation. 

2. Turn rate was approximately the same, i.e., .007 rads/ 
sec for both situations. However, as might be expected from 
a human, the turn rate exhibits a much more sinusoidal be- 
havior as compared to the optimal solution which immediately 
sought and settled down to a steady state value. 

3. The state variables under real time simulation showed 
a more pronounced fluctuation. 

4, An advantage of the real time simulation over the 
optimal solution was that the human operator in all cases, 
1.e€. both thrust failure and thrust reversal, was able to 
maneuver the ship back to the original track, whereas the 
optimal solution maneuvered the ship to a parallel track and 
never returned to the original. This, of course, stems from 
the method used to find the "steady state" values around 
which the equations of the system were linearized. 

Nevertheless, the difference of 18 feet between the original 
track and the one that is reached at steady state by the opti- 
mal algorithm is not significant compared either to the ship's 


wen (100 feet) or to the tactical considerations of the ship. 
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mm NONLINEAR VERSUS LINEAR SYSTEM RESULTS 

In order to determine how "optimal" the response of the 
nonlinear system was with the use of the state feedback 
values obtained from the linear regulator design, a com- 
parison of values of the states and the performance indices 
were made during a failure. 

The results are shown in Table ITII. As can be easily 
verified, the differences in steady state values obtained 
from the linearized system and the nonlinear system are 


ca) 


Maite small, ranging from .77% to 1.77%, with the exception 
of the steady state thruster angle (15% difference). 
The above results point out that the linear assumptions 


and approximations used in this study hold quite well over 


the operating range. 


eB IE iy 
Comparison of Linear and Nonlinear Models 


LINEAR NONLINEAR 


Surge velocity (ft/s) i Se Jie 74.63 
Sway velocity (ft/s) 1.88 Ley, 
Angular velocity (rad/s) wales ~.211x107° 
Yaw angle (rad) ~,0248 -.025 
Pusplacement in x direction (ft) 7453 7342 
Displacement in y direction (ft) HGreeo7 0.88 
Steady state thruster angle (rad) .15256 .1804 


81 





Vi. CONCLUSIONS AND RECOMMENDATIONS 


me CONCLUSIONS 

Comparison of the results of all three methods examined 
leads to the conclusion that the optimal control solution is 
Dy far the best one. 

Control effort, be it thruster angle or thruster magni- 
tude variation, is less with the optimal control solution. 
his is of great consideration not only from the aspect of 
wear and tear of the mechanical parts but also from the 
point of ship performance as a naval unit, especially since 
its speed is so far greater than the conventional ships. 

Augmentation of human operator control inputs with the 
optimal solution during failure scenarios presents the best 
solution in order to maintain positive control of the ship 
under any operating conditions. This approach would combine 
the best attributes of the human operator and the optimal 


control schemes. 


B. RECOMMENDATIONS 

Based on the previous analysis, the following could be 
investigated in the future. 

1. A sufficiently complex constrained parameter optimi- 
zation routine (COPES-CONMIN) can be used to arrive logically 
at a set of relative system response with respect to physical 
Constraints. This optimization will still be dependent on the 


performance index selected by the design engineer and his 
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Mrs ability to calculate explicitly its value for a given 
dynamic response [Ref. 11]. 

2. The right-hand side of the equations of motion could 
be expanded to include more terms to provide an increase in 
quantitative accuracy of the drag terms such as bubble drag 
and include also a constant power propulsion plant versus a 
constant thrust one. This could be done easily since the 
terms are additive. 

3. Development of a computer program to augment the 
control efforts commanded by a human operator should be 
investigated, tested and implemented on the RTSSD of Ref. 4. 

u, It might prove worthwhile to investigate more closely 
the trial and error design for improved response character- 


maeress or this method of control. 





APPENDIX A 


LINEAR REGULATOR DESIGN THEORY 


Optimal control theory provides a means for utilization 
of state variable feedback in order to control a system by 
minimizing a performance index or cost function [Refs. 7 and 


8]. 


Having a linear system defined as x = Ax + Bu where u 


~e ~e 


1s the control input vector, the linear op imal control law 


[m= detined as opt? = -Gx(t). 


The assumptions under which this solution is valid are: 
meine eomenrol vector as uneonstrained. 


2. The performance index is defined as 


= je {x" Qx + uw Ry} at 
oO 


where Q is a positive semidefinite or positive definite real 
Symmetric matrix and R is a positive definite real symmetric 
Matrix. The matrices Q and R determine the relative impor- 
tance of the square of the error and energy expenditure 
during the transient period. 


3. The system is controllable and state observable. 


With these assumptions, substitution for uw leads to 


4 

ey (x Qx + G(t)'x R G(t)x) dt 
O 

J = ftt{x’ (Q + G(t)” R G(t))x} at 


@) 
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The matrix G(t) is found to be 


digi 


eek BE P(t) 


~ 


Here P(t) satisfies the matrix Riccati equation 


=) 


P(t) + Q - P(t) BRB P(t) + P(t)A + A'P(t) = 0 


subject to the boundary condition P(tf) =i 

The elements of the matrix G(t) are referred to as the 
feedback coefficients, since the optimal control consists of 
a time-weighted linear combination of the state variables. 

This result is important since it indicates that in an 
optimal system all the state variables are to be fed back, 
not only the output, as it is customary in classical control 
theory. A graphical representation of this optimal system 
is presented in Fig. 42. 

The Riccati equation is integrated backwards in time from 
the known terminal condition P('~) = 0 until the desired 
degree of accuracy 1s reached. When the final time approaches 
Marainity, the P Matrix becomes constant, thereby simplifying 
its determination. 

As can be seen for the previous presentation, a different 
linear optimal control is generated for every relative set 
of weighting matrices, Q and R, provided. Additionally, the 
system response is optimal in a mathematical sense only, since 
no assurance on acceptable system overshoots and input peak 
Values can be guaranteed. 


A minimum area could have an excessive peak which would 
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| 
provide a physically unrealizable control scheme. The 
design synthesis thus requires an initial guess at the rela- 
| tive importance multipliers, followed by subsequent adjust-- 


ments based on time simulation of the controlled system 


dynamics . 
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